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(54) Phase Estimation And Compensation .n Orthogonal Frequency Division Multiplex (OFDM) 
Systems 

(57) Methods are presented for phase estimation 
and compensation in Orthogonal Frequency Division 
Multiplex systems. The methods include the step of ex- 
tracting training tones from a received digitized data 
burst. Some of the methods further include the steps of 
determining a channel impulse response based on the 
extracted tones; estimating phase shift caused by com- 
mon phase noise; correcting the impulse response and 
performing a Fast Fourier transform to provide a chan- 
nel estimate. Some of the methods include estimating 
the phase shift caused by common phase noise before 
an average impulse response is updated with the cur- 
rent impulse impulse. Some of the methods include fre- 
quency shifting the spectrum of the channel Impulse re- 
sponse and providing partial estimates. 
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Description 

FIELD OF INVENTION 



5 [0001] This invention relates to Orthogonal Frequency Division Multiplex (OFDM) systems and more particularly to 
phase estimation and compensation. 

Background of Invention 

w [0002] Orthogonal Frequency Division Multiplex (OFDM) is used for the air interface for the IEEE 802.11a, ETSI 
Hiperlan-ll, terrestrial digital video broadcast (DVB-T), and Broadband Wireless Internet Forum (BWIF) standards. It 
is also used in wireline transmission, notably Asynchronous Digital Subscriber Lines (ADSL). The OFDM technique 
sends many carriers (sometimes thousands of carriers) in parallel on adjacent frequencies within a frequency band. 
The frequencies are variously called frequency "bins", tones, or subbands. Tones is the term used in the following 

is description 

[0003] The OKDM technique relies on the orthogonality properties of the fast Fourier transform (FFT) and the inverse 
fast Fourier transform (IFFT) to eliminate interference between the carriers. The FFT can be thought of as a bank of 
filters, with each filter having the identical frequency response as the other filters, but centered on different center 
frequencies. The response of each filter in the filter bank is a maximum at the center frequency of the filter, but zero 
20 at the center frequencies of the other filters. Thus, if the carriers are centered on the center frequencies of the filters, 
intercarrier interference is eliminated. If a carrier frequency is not perfectly centered, then its frequency is not at the 
zero response points of the FFT filters, and some of the carrier energy appears in the FFT outputs at the other fre- 
quencies. In FFT nomenclature, this is called leakage. Leakage is the main source of intercarrier interference in OFDM 
systems. 

25 [0004] At the transmitter, the precise setting of the carrier frequencies is performed by an I FFT. The data is encoded 
into constellation points by multiple (one for each carrier) constellation encoders. The complex values of the constel- 
lation encoder outputs are the inputs to the IFFT. For wireless transmission, the outputs of the IFFT are converted to 
an analog waveform, unconverted to a radio frequency, amplified, and transmitted. 

[0005] At the receiver, the reverse process is performed. The received signal is amplified, downconverted to a band 
30 suitable for analog to digital conversion, digitized, and ultimately processed by a FFT to recover the carriers. The 
multiple carriers are then demodulated in multiple constellation decoders (one for each carrier), recovering the original 
data. 

[0006] Since an I FFT is used to combine the carriers at the transmitter and a corresponding FFT is used to separate 
the carriers at the receiver, the process has potentially zero intercarrier interference (ICI). However, practicalities such 
35 as synchronization errors, frequency tracking errors, and phase noise introduce intercarrier interference by shifting or 
jittering the carriers from the center frequencies of the filters and allowing some leakage. Some of these effects can 
be compensated, and some cannot. This patent is concerned with compensation methods for phase noise degrada- 

[0007] In wireless transmissions, the signal can be reflected or scattered from buildings, vehicles, trees, vegetation, 
40 and terrain features. Multiple copies of the signal, each with a different time delay that depends upon the path traversed, 
are summed at the receive antenna, This phenomenon is called multipath transmission, or more simply multipath. 
Multipath causes fading and attenuation in the frequency band, which if uncompensated would cause unacceptably 
large numbers of errors in the decoding process. 

[0008] To avoid this, carriers with known amplitude and phase are transmitted for the purpose of measuring the 
45 wireless transmission channel. These earners are known as training tones or pilot tones. In some systems, they are 
transmitted in a separate burst composed entirely of training tones. In other systems, such as the BWIF standardized 
system to be described herein, they are interspersed among the data carriers (data tones.) In either case, since the 
training tones are known apriori, the response of the wireless channel at the training tone frequencies is easily deter- 
mined. Once this is done, the channel responses at the data tone frequencies can be interpolated from the known 
so responses at the training tone frequencies. This process is called channel estimation. The measured and interpolated 
channel responses, known as channel estimates, are used in the decoding process to minimize decoding errors. 
[0009] Phase noise is introduced into the signals by phase noise in the local oscillators used to upconvert and down- 
convert the signal, as well as by timing jitter of the sampling in the digital to analog and analog to digital converters at 
the transmitter and receiver, respectively. Each carrier, which would have a line spectrum otherwise, now has a phase 
55 noise spectrum associated with it. 

[0010] For OFDM systems, phase noise can be classified into two categories - common phase noise and foreign 
phase noise. Common phase noise is so named because it affects all carriers equally. It is generated by the local 
oscillators 1 phase noises and the data converters' timing errors, which modulate each carrier equally. This causes a 
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phase rotation on each and every carrier that is equal in sign and amount. Foreign phase noise, on the other hand, is 
caused by the leakage of the phase noise spectra of the other carriers ("foreign" carriers) into the filter response of a 
given carrier. Foreign phase noise appears as random noise from carrier to carrier. The effects of phase noise on 
OFDM systems are analyzed in the article "Influence of RF Oscillators on an OFDM Signal", by Claus Muschallik. This 

5 was published in the IEEE Transactions on Consumer Electronics, Vol. 41 , No. 3, August 1995, pp. 592-602. 

[0011] Fig. 1 illustrates a functional block diagram of a wireless receiver system 10 employing a channel estimator 
20. A data signal or burst is received by an antenna 14, which picks up the data signal and couples the data signal to 
front end processing 12 . The data signal or burst can include training and calibration information. The processing 12 
amplifies the data signal, converts the data signal to an IF frequency and filters the data signal to remove signal outside 

10 the desired frequency band. The front end processing 12 is coupled to one or more analog-to-digital converters 16 
that sample the data signal and provide digitized signal output. The digitized signals of the one or more analog to digital 
converters are provided to a digital preprocessor 18. The digital preprocessor then performs a FFT on the digitized 
signal. The FFT on the digitized signal converts the signal from the time domain to the frequency domain so that the 
frequencies or tones carrying the data can be provided. The frequencies and tones can then be demodulated or de- 

15 coded. The demodulation of the tones requires information relating to the wireless channel magnitude and phase at 
each tone. The effects of the dispersion caused by the channel needs to be compensated prior to decoding of the 
signal, so that decoding errors can be minimized. A channel estimator 20 Is provided to determine the amount of phase 
rotation and magnitude perturbation applied to the tones by the channel. Since the training tones are transmitted with 
known magnitude and phase, the channel response at the training tones is easily determined. The known channel 

20 response at the training tones is then interpolated in the frequency domain to determine the channel response at the 
data tones. A cyclic interpolation procedure can be employed. The channel estimate is then provided to the data de- 
modulator 22 for demodulation of the digital data signal. The demodulated data signal is then transmitted to a data 
postprocessing block 24 for further signal processing. The data post processing performs error correction utilizing the 
information from the data demodulator 22. The channel estimator 20 extracts the training tones from the number of 

25 tones in the data signal or data burst and performs several signal processing steps on the training tones. After correcting 
for any transmitted magnitude and phase differences, the channel estimator 20 performs an IFFT on the training tones 
to provide a channel impulse response. The channel impulse response will eventually be zero padded and have a FFT 
performed on it to obtain channel estimate values at frequencies between the training tone frequencies. 
■ s [0012J The BWIF channel estimation method is shown in Figure 2. The procedure starts by providing a set of fre- 

30 quency samples of training tones H(k) and at block 31 extracting the training (also called pilot) tones Hj^k) from the 
'■ set of all tones H. Note that N denotes the number of all tones, and N-n- denotes the number of training tones. Modulation 
on the training tones is removed at block 32, and the edge tone, which is transmitted at reduced amplitude, is re-scaled 
to proper amplitude at step 33. In block 34, an inverse fast Fourier transform (IFFT) on the training tones is performed 
to obtain the channel impulse response. The impulse response is scaled to unity in block 35. In block 36 the impulse 

35 response is averaged with previous bursts to reduce noise. The average impulse response is then zero padded (block 
37) and an N-point FFT is performed at block 38 to produce the final channel estimate. Note that it is possible to omit 
the averager at block 36 (or equivalently, set a equal to zero in the averaging equation in block 36). 
[001 3] Common phase noise causes ail of the tones (H and H-py) to be phase sh ifted by a common amount. Therefore, 
every term of the impulse response is also phase shifted by the same amount. This phase shift is a zero mean random 

40 variable and changes for each burst. The problem occurs at block 36 of Figure 1 , when multiple impulse responses 
are averaged to produce an average impulse response. The averaging is intended to reduce the noise in the final 
estimate. However, it also averages the common phase noise to near zero. This can be a problem in that the data 
tones will have the common phase shift, but the channel estimate produced from the average impulse response will 
not. This resurts in increased error rates in the demodulation of the data. For the standard BWIF processing, this is 

45 only a problem if the averager is enabled (a *0 ). If the averager is not enabled, the common phase noise is present 
identically in both the channel estimate and the data, and is thus compensated in the standard demodulation process. 
[0014] In addition to this problem, Texas Instruments has developed (see Appendix A; TI-32977, Serial No. 
10/001 ,986 filed October 31 , 2001) a computationally more efficient method of obtaining channel estimates for the 
BWIF system. This application is incorporated herein by reference. The Tl method estimates 1/K of the tones (typically 

so 1/4 of the tones) at each burst, and uses K bursts (typically 4 bursts) to estimate the tones. At each burst, a channel 
estimate for ail tones is obtained by using the N/K tone estimates computed on the current burst, and the (K-1) # N/K 
tone estimates from the (K-1 ) previous bursts. 

[001 5] Phase noise causes two problems here. First, any impulse response averaging will average the phase noise 
to zero, just as in the original BWIF method. Since the data will be phase shifted (by the common phase noise), but 
55 the channel estimate is not, there is an error introduced in the demodulation of the data, which causes increased error 
rates. Secondly, since the complete channel estimate is formed over K bursts, each N/K set of estimates was obtained 
at a different phase state (from the common phase noise). This makes the channel estimates obtained with this method 
noisy. 
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SUMMARY OF INVENTION 

[001 6] In accordance with one embodiment of the present invention a metnod for estimating the burst-to-burst phase 
shift and adjusting the channel estimates to compensate for It is provided by first obtaining an estimate of the phase 
5 shift caused by the common phase noise and second by either phase correcting the impulse noise or tones in the 
channel response. 

[0017] The actual implementation procedure of these steps is specific to the channel estimation procedure. One 
embodiment is BWIF standard channel estimation with impulse response averaging. 
[0018] A second is a Texas Instruments Inc. estimation procedure with impulse response averaging. 
10 [0019] A third is the Texas Instruments Inc. channel estimation without impulse response averaging. 

[0020] A fourth embodiment is an alternative implementation for the BWIF standard to estimate the phase angle 
based on the value of h avg before it is updated with the current impulse response. 

[0021] A fifth embodiment is an alternative embodiment using the Tl estimation method with impulse response av- 
eraging incorporating the phase estimation method of the fourth embodiment to estimate the phase angle based on 
is the value of h avg before it is updated with the current impulse response. 

DESCRIPTION OF DRAWINGS 

[0022] Figure 1 illustrates a functional block diagram of a wireless receiver employing a channel estimator. 
20 [0023] Figure 2 illustrates BWIF channel estimation procedure according to the prior art. 

[0024] Figure 3 illustrates BWIF channel estimation procedure incorporating phase noise estimation and compen- 
sation according to one embodiment of the present invention. 

[0025] Figure 4 illustrates an embodiment of the Texas Instruments Channel Estimation Procedure with impulse 
response averaging. 

25 [0026] Figure 5 illustrates the Texas instruments estimation method with impulse response averaging incorporating 
phase estimation and compensation. 

rO027] Figure 6 illustrates the Texas Instruments Channel estimation routine without impulse response averaging. 
[0028] Figure 7 illustrates the Texas Instruments Channel Estimation procedure ( without impulse response averag- 
ing) incorporating phase estimation and compensation. ... 
30 [0029] Figure 8 illustrates the standard BWIF channel estimation procedure incorporating an alternate implementa- 
tion of phase noise estimation and compensation. 

[0030] Figure 9 illustrates the Texas Instruments channel estimation method with impulse response averaging incor- 
porating the alternate phase estimation and compensation. 

35 DESCRIPTION OF PREFERRED EMBODIMENT OF THE PRESENT INVENTION 

[0031] This patent addresses means to compensate for common phase noise in BWIF standardized and similar 
OFDM systems. These systems have training tones interspersed among the data tones on each burst. There are two 
broad categories to consider. In these systems, the channel estimates can be obtained either from the training tones 

40 of the current burst exclusively (Case 1 ), or from an average of the training tones in the current and previous bursts 
(Case 2) On any given burst, common phase noise has affected the tones by introducing a phase rotation, common 
to all the tones. In Case 1, since the channel estimates are obtained from these phase shifted training tones, the 
channel estimates are similarly phase shifted, and the phase shift on the data is automatically compensated by the 
phase shift on the channel estimates in the decoding process. 

45 [0032] In Case 2, If the channel estimates are obtained from an average of training tones over several bursts, then 
the averaging process has likely reduced the phase rotation on the channel estimates. The data tones however have 
the phase shift of the current burst. Since the channel estimates and the data do not match in phase shift, some beneW 
can be obtained by a compensation procedure. In the procedure to be described herein, the difference in phase shift 
between channel estimates and the current burst is first estimated. Then, the estimate is used to rotate the channel 

so estimate to the phase state of the current burst prior to the data decoding operation. 

r0033] In addition to these cases, Texas Instruments (Tl) has developed several channel estimation procedures that 
are computationally more economicalthanthestandardBWIFestimalionprocedure. These are described .n ^plication 
serial no 10/001 986 filed October 31 of 2001 (TI-32977) entitled "A Computationally Efficient System and Method for 
Channel Estimation." This application is incorporated herein by reference. A copy is provided in Appendix A. These 

55 economical procedures obtain channel estimates for different tones on successive bursts. Thus, the phase rotation 
introduced into the channel estimates by the common phase noise is different for the different tones, depending upon 
which buistthe estimate was obtained. Phasenoisecompensationproceduresforthese procedures are also described 
herein. 
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[0034] The solution is a two step process. First, we obtain an estimate of the difference between the common phase 
noise induced phase shifts of channel estimates and the phase shifts of the current burst. Second, we phase align 
either the impulse response or the tones In the channel estimate with the phase of the current burst. The actual imple- 
mentation of these steps is specific to the channel estimation procedure. There are five embodiments. 

5 

1 . The BWIF standard channel estimation method with impulse response averaging. 

2. The Tl channel estimation with impulse response averaging. 

3. The Tl channel estimation without impulse response averaging. These will each be discussed below. 

4. An alternative implementation for the BWIF standard including an alternative implementation of phase noise 
10 estimation and compensation. 

5. The Tl channel estimation method with impulse response averaging incorporating the alternate phase estimation 
and compensation. 

1 , Phase noise compensation for the BWIF standard channel estimation procedure with impulse response averaging 

15 

[0035] Figure 3 illustrates the first embodiment ot the present invention. The procedure starts as done previously in 
Figure 2 by providing a set of frequency samples of training tones H(k) and at block 41 extracting the training (also 
called pilot) tones Hj^k) from the set of all tones H. The modulation on the training tones is removed at block 42, and 
the edge tone, which is transmitted at reduced amplitude, is re-scaled to proper amplitude at step 43. In block 44, an 
20 inverse fast Fourier transform (IFFT) on the training tones is performed to obtain the channel impulse response. The 
impulse response is scaled to unity in block 45. In block 46 the impulse response is averaged with previous bursts to 
reduce noise. 

[0036] The modifications to the BWIF standard channel estimation procedure of Figure 2 are shown in Figure 3 in 
bold outline. Here, all of the current impulse response samples h(n) have been phase shifted due to the phase noise, 
25 but the average impulse response h avg (n) has near zero phase shift. To estimate the difference in phase between the 
current impulse response and the average impulse response, the following computation is performed in block 47: 



30 
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where O 0 is the estimate of the phase difference between the average impulse response and the current impulse 
35 response. Since the data is phase-aligned with the current impulse response, but the channel estimate is obtained 
from the average, the average impulse response is then phase corrected by - «> 0 to phase-align it with the current data. 
This is accomplished with the following computation: 



A'avg(n) = ha Vg (n).e ^° » n= 0,1,— , Ntt - 1 

(2) 



[0037] Here, h' is the phase corrected average. The average impulse response is then zero padded (block 48) 
using h' avg and an N-point FFT is performed at block 49 to produce the final channel estimate. 

so 2. Phase noise compensation for the Tl channel estimation procedure with impulse response averaging. 

[0038] The Tl channel estimation procedure with impulse response averaging is shown in Figure 4. Blocks 51 through 
56 are identical to those in the standard BWIF channel estimation method (Blocks 31 -36) shown in Figure 2. However, 
in the remaining blocks, the processing is altered so that only 1/K of the tones are estimated at any one burst. The 
55 channel estimate for a burst is composed of the N/K tones estimated currently and (K-1 )N/K tones estimated at previous 
bursts Block 57 is a burst counter, modulo K. The value of b is a burst count index used for later computations. At 
block 58, the average impulse response is frequency shifted or "twiddled" by multiplication of h avg (n) wi th enW. 
At block 59, the twiddled average impulse response is zero padded to length N/K. At block 60, an N/K point FFT is 
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performed to obtain the channel estimate for the N/K tones being estimated on the current burst. These are combined 
with the remaining tones, estimated on the previous bursts, to form the complete (all tones) channel estimate at block 
61 In block 61 , the notation H.., refers to the channel estimate at the previous burst. A more detailed descnption is 
found in Appendix A of the present application. This is described in Application serial no. 10/001 .986 (TI-32977) filed 
October, 31 of 2001 entitled "A Computationally Efficient System and Method for Channel Estimation. This application 
is incorporated herein by reference. 

[0039] The impulse response averager at block 56 causes the same problem here that it caused in the standard 
BWIF channel estimation procedure. That is, it averages the common phase noise to near zero. Figure 5 shows the 
modifications to Figure 4 required to estimate and compensate for the phase shift induced by common phase noise. 

10 The blocks 62-67 remain the same as in Figure 4. Block 68 is added to the output of block 67. Blocks 69 - 72 are like 
blocks 57-60 in Figure 4. Block 73 is modified from block 61 in Figure 4. The blocks that are added to or modified from 
Figure 4 are in bold outline. The phase estimation and correction block 68 in Figure 5 estimates the phase difference 
between the average impulse response and the current burst, then applies a phase correction to the average to bring 
it "into phase" with the currentburst. The modified average is then used in block 70 (instead of the unmodified average). 

is The result is that the channel estimates computed at block 71 are properly phased with respect to the common phase 
noise of the burst and the data. 

,00401 Further modifications are required at block 73, where the current bursts estimates are combined with previous 
burst estimates to obtain the complete (all tones) estimate for the burst. The tones of the previous burst were phase 
aliened wilh the previous burst, not the current burst Thus, a correction term on the previous burst's estimates is 
so required to match the phase of the previous burst with that of the current burst. This is done by multiplying the previous 
burst tones with ,-<+ -* •> . Here, and O., are the phase estimates at the current burst and previous burst, respectrvely. 

3. Phase noise compensation for the Tl channel estimation .procedure when impulse response averaging is not used 

25 r0041l There may be circumstances where it is advantageous to disable the impulse response averaging and com- 
pute channel estimates based on the impulse responses of a single burst. The poses a different scenario for phase 
estimation and compensation. . 
r00421 Figure 6 shows the Tl channel estimation procedure without impulse response averaging. Steps 81 through 
85 are ! idenfical to those of the BW.Fstandard procedure by blocks 31 -35 shown in Figure 2. However, '"^defusing 

30 the impulse response of step 85 to compute the channel estimates for all tones, the Tl procedure uses the impulse 
response to compute the estimates for a fraction of the tones (N/K tones, where K is an integer power of 2). Since the 
impulse response is not averaged wrth previous impulse responses, it carries the phase shifts caused by < ommon 
ohase noise Block 86 of the procedure increments a burst counter (modulo K) as prev.ously. The index b is used in 
tcompu^ 

response is multiplied by twiddle factors a*",*. The result is zero padded to length WK n block 8 block 89 ^ 
N/K point FFT is performed to obtain the channel estimates at the N/K frequences. This is combined with estimates 
of the other frequencies (estimated at previous bursts) in block 90. „„„„„„ 
[0043J The channel estimates obtained at block 89 are all uniformly shrfted by the phase shift mduced by common 
Phase noise. However, the tones estimated at previous bursts and used in block 90 have a different phase shift. So, 
the estimate obtained at the output of block 90 will have an incorrect phase shift on (K-1 N/K ^ ue " cies - 
00441 To correct this, the processing is altered as shown in Figure 7. The processing in blocks 91-95 » like tha in 
blocks 81-85 in Figure 6. Blocks 97-100 are like blocks 86-89 in Figure 6. At block 96 a phase estimation eshmates 
the difference in phase between the present burst and the previous burst. This is computed as 



35 



40 



43 



SO 




55 



where h , (n) are the samples of the previous burst's impulse response. 
[00451 TheangleV-stheestimatedphased^^ 

estimate H., , where the subscript (-1 ) indicates the previous burst. Thus, all the tones in the total estimate H have the 
same phase with respect to phase noise. 
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4. An alternative implementation for the BWIF standard including an alternative implementation of phase noise 
estimation and compensation. 

[0046] Referring to Figures 8 and 9 there is an alternative procedure for compensating for phase noise. The alter- 
5 native implementation is shown in Figure 8forthe BWIF standard processing and is shown in Figure 9 fortheTI channel 
estimation procedure with impulse response averaging. Blocks 111-115 and 118 and 119 in Figure 8 is the same as 
the standard BWIF blocks 31-35 and blocks 36-38 in Figure 2 The Idea behind the alternative phase compensation 
procedure is to estimate the phase angle & 0 based on the value of h avg (block 1 1 6)before it is updated with the current 
impulse response. Then, when the average is updated, the old average is phase corrected by multiplication by c -A to 
10 "phase align" it with the current impulse response. This is shown in block 11 7 of Figure 7. With this alternative approach, 
the average impulse response is phase aligned with current impulse response measurement. 

5. The Tl channel estimation method with impulse response averaging incorporating the alternate phase estimation 
and compensation 

15 

[0047] Figure 9 illustrates a similar modification to the Tl channel estimation method as described in Figure 4. Blocks 
121 -125 correspond to blocks 51-55 in Figure 4 and blocks 128-132 correspond to blocks 57-61 in Figure 4. Block 126 
estimates the phase angle <I> 0 based on the value of h avg . The old average is phase corrected by multiplication by c"* 
to "phase align" it with the current impulse response. This is shown in block 127 of Figure 9. 
20 [0048] Although preferred embodiments have been described, it will be apparent to those skilled in the relevant art 
that various modifications, additions, substitutions and the like can be made without departing from the spirit and of 
the invention and these are therefore considered to be within the scope of the invention as defined in the following 
claims. 

25 

Claims 

1 . A method of channel estimation comprising the steps of: 

30 extracting training tones from a received digitized data burst; 

determining a channel impulse response based on the extracted tones; 
: estimating phase shift caused by common phase noise; 

correcting the impulse response; and 

performing a Fast Fourier Transform to provide a channel estimate. 

35 

2. The method of Claim 1 including the step of zero padding the impulse response or tones and performing a Fast 
Fourier Transform on the zero padded average impulse response to provide a channel estimate. 

3. A method of channel estimation comprising the steps of: 

40 

extracting training tones from a received digitized data burst; 
determining a channel impulse response based on the extracted tones; 
averaging the impulse response with previous impulse responses; 

estimating the difference in phase between the current impulse response and the average impulse response; 
45 correcting the average impulse response to phase align it with current impulse response; 

and 

performing a Fast Fourier Transform on the average impulse response to provide a channel estimate. 

4. The method of Claim 3 including the step of zero padding the impulse response or tones and performing a Fast 
so Fourier Transf orm on the zero padded average impulse response to provide a channel estimate. 

5. A method of channel estimation comprising the steps of: 

extracting training tones from a received digitized data burst; 
55 determining a channel Impulse response based on the extracted tones; 

averaging the impulse response with previous impulse responses; 

estimating the difference in phase between the current impulse response and the average impulse response; 
correcting the average impulse response to phase align it with current impulse response; 
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frequency shifting the spectrum of the channel Impulse response; 

performing a Fast Fourier Transform on the frequency shifted channel impulse response to provide a partial 
channel estimate; and 

combining the partial channel estimate with at least one other partialchannel estimate to form an aggregated 
channel estimate. 

6 The method of Claim 5 including the step of zero padding the Impulse response and performing a Fast Fourier 
Transform on the zero padded frequency shifted channel impulse response to provide a partial channel estimate. 

7. A method of channel estimation comprising the steps of: 

extracting training tones from a received digitized data burst; 
determining a channel impulse response based on the extracted tones; 
averaging the Impulse response with previous impulse responses; 

estimating the difference in phase between the current impulse response and the average impulse response; 
correcting the average impulse response to phase align it with current impulse response; 
frequency shifting the spectrum of the channel Impulse response; 

performing a Fast Fourier Transform on the frequency shifted channel Impulse response to provide a partial 
channel estimate; 

phase correcting at least one previous partial channel estimate; and 

combining the partial channel estimate with at least one other previous phase corrected partial channel esti- 
mate to form an aggregated channel estimate. 

8 The method of Claim 7 including the step of zero padding the impulse response and performing a Fast Fourier 
" Transform on the zero padded frequency shifted channel impulse response to provide a partial channel estimate. 

9. A method of channel estimation comprising the steps of: 

extracting training tones from a received digitized data burst; 
determining a channel impulse response based on the extracted tones; 
estimating the difference in phase between the current burst and previous burst; 
frequency shifting the spectrum of the channel impulse response; 

performing a Fast Fourier Transform on the frequency shifted channel impulse response to provide a partial 

comb^ estimate with at least one other partial channel estimate to form an aggregated 

channel estimate. 

1 0 The method of Claim 9 including the step of zero padding before performing a Fast Fourier Transform on the zero 
padded frequency shifted channel impulse response to provide a partial channel estimate. 

11. A receiver or system implementing the method of any preceding claim. 
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